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We report first-principles all-electron density-functional based studies of the electronic structure, 
magnetic ordering and anisotropy for the V15 molecular magnet. From these calculations, we de- 
termine a Heisenberg Hamiltonian with five antiferromagnetic and one ferromagnetic exchange cou- 
plings. We perform direct diagonalization to determine the temperature dependence of the suscep- 
tibility. This Hamiltonian reproduces the experimentally observed spin 5=1/2 ground state and 
low-lying 5=3/2 excited state. A small anisotropy term is necessary to account for the temperature 
independent part of the magnetization curve. 
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With the continued interest in the fabrication and op- 
timization of miniaturized magnetic devices future 
design considerations will require an understanding of 
nanoscale magnetic systems. In order to transition such 
materials into simple devices it is necessary to be able 
to explain how interactions such as spin-spatial coupling 
and spin-spin exchange effects may couple collectively to 
create a seemingly single-spin system. Further, it is nec- 
essary to determine the temperature range at which such 
systems will indeed behave collectively. In general the 
properties of a nanoscale system of coupled spins depend 
directly on the strength of the exchange-parameters and 
on the size and sign of the anisotropy energy due to spin- 
orbit coupling. While these parameters are generally de- 
termined by the transition metal atoms, the ligands and 
other nonmagnetic atoms are responsible for stabilizing 
the array of spins. Requisite to a complete computa- 
tional understanding of such spin systems is the ability 
to account for the strong ligand-metal interactions and 
to determine whether the behavior of a given spin sys- 
tem is mainly mediated by the anisotropy, by spin-spin 
coupling or by a combination of the two. 

Recently, the Mni2-Acetate and Feg molecules have 
attracted considerable interest because they behave as 
high-single-spin systems (total spin 5=10) at tempera- 
ture ranges on the order of 20-60K. Due to their large 
magneto-molecular anisotropy energy these systems re- 
tain their moment orientation at reasonably high temper- 
atures and exhibit the phenomena of resonant tunneling 
of magnetization at well defined magnetic fields . 

The K6[Vi5As6042(H20)] 8H2O molecular crystal, 
first synthesized by Miiller and Doring |^,^, represents 
a transition-metal spin system in the same size regime 
as the Mni2 and Feg molecular crystals. In contrast to 
Mni2 and Feg molecules, the V15 molecule is thought 
to behave as a weakly anisotropic magnet composed of 
15 spin s=l/2 particles which couple together to form a 
molecule with a total spin 5=1/2 ground state. Besides 
the fundamental interest in understanding quantum ef- 



fects in these nanomagnets they might be also relevant 
for implementations of quantum computers 0. Calcu- 
lations on such correlated systems present a challenge 
to mean-field frameworks such as density-functional the- 
ory because it is often not possible to construct a sin- 
gle collinear reference state which preserves the inherent 
symmetry of the system and has the correct spin quan- 
tum numbers. 

This work utilizes an efficient coupled multilevel anal- 
ysis which relies on fitting density-functional energies to 
mean-field Heisenberg or Ising energies in order to de- 
termine the exchange parameters. The approximate ex- 
change parameters gleaned from the first N Ising config- 
urations were used to find the next lowest energy Ising 
configuration and subsequently to improve the param- 
eterization of the exchange parameters. "Self Consis- 
tency" of this approach is determined when the predicted 
Ising levels are unchanged by the addition of data from 
new Ising configurations. The coupling of the density- 
functional method to a classical Ising representation al- 
lowed us to determine the exchange parameters by con- 
sidering only several spin configurations. Further, the 
resulting ground-state spin configuration within density- 
functional theory exhibits the correct spin projection of 
1/2. With the exchange parameters determined, we di- 
agonalize the complete many-body Heisenberg Hamilto- 
nian to calculate the susceptibility and spin correlation 
functions for comparison with experiment. The many- 
body basis is complete, so all states are allowed including 
non-coUinear spin arrangements and quantum disordered 
phases js). 

Starting from X-ray data we generated several unit 
cells and isolated a single K6[Vi5As6042(H20)] unit. In 
order to optimize the geometry within the quasi- 1?3 sym- 
metry of the Vi5 molecule we initially replaced the 
statistically oriented water molecule at the center of the 
molecule by a neon atom and used an 5=3/2 spin config- 
uration which does not break the crystallographic sym- 
metry. The Ising configuration of this molecule consists 
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of three aligned spin-i V atoms in the triangle and equiv- 
alent upper and lower hexagons composed of a ring of 
antiferromagnetically (AF) coupled spin-i V atoms. 

The geometry of the molecule was then optimized 
within the all-electron density-functional methodology 
using the generalized-gradient approximation (GGA) 
pO| . The calculations were performed with the Naval 
Research Laboratory Molecular Orbital Library (NRL- 
MOL) Calculations on 49 geometrical configura- 

tions were performed during the conjugate-gradient re- 
laxation of the molecule. Subsequent calculations show 
that the geometrical, electronic and magnetic properties 
of this system are unaffected by the presence or type of in- 
ert moiety enclosed within the void. Using this geometry, 
we performed eleven additional calculations on different 
spin configurations (See Tabic ||) to determine the six ex- 
change parameters ( J's) of the Hciscnbcrg Hamiltonian 

H = ^JijSi-Sj, (1) 

as well as the spin configuration of the density-functional 
ground state The J's used in the above Hamiltonian 
are defined according to Fig. |l|. As shown in Table |[ we 
have included high-spin configurations (XI, XII, XIII), 
which generally have some symmetry as well as lower- 
spin non-symmetric configurations. The energy for the 
high-spin 5 = 15/2 ferromagnetic (FM) state (XIII) of 
873 meV is predominantly caused by a large AF exchange 
coupling (J) between the most closely bonded hexagonal 
V atoms. However, the 113 meV splitting between the 
5=9/2 and 5=15/2 states (XII and XIII) shows that 
there is a reasonably strong AF coupling, approximately 
18 meV on average, between the triangular and hexago- 
nal atoms. All of the data displayed in Table | has been 
used to determine the exchange parameters from a least 
square fit to the mean-field solution of the Heisenberg 
Hamiltonian (^ . The fit is very good with errors ranging 
from 0.1-1.55 mcV. The fit leads to exchange parameters 
of J = 290.3 meV, J'= -22.7 meV, J"=15.9 meV, Jl = 
13.8 meV, J2 = 23.4 meV and J3 = 0.55 meV, where pos- 
itive numbers correspond to AF and negative ones to FM 
interactions. The ferromagnetic interaction J' is a sur- 
prising result and deserves further discussion since it 
is qualitatively different from earlier assumptions based 
on entirely AF interactions |^,^. Ferromagnetic cou- 
pling is possible without polarizing the oxygens through 
a 4'th order process similar to super-exchange. In super- 
exchange, the intermediate state has the lowest c?-orbital 
on the V doubly occupied with up and down electrons 
. However, electrons can also hop to higher energy d- 
orbitals on the Vs. In this case both parallel and antipar- 
allel spins are allowed without violating the Pauli exclu- 
sion principle, but Hund's rule coupling prefers parallel 
alignment. The super-exchange process (same c?-orbital) 
completely excludes the process with same-spin electrons 
while the FM process (different d-orbitals) merely favors 
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FIG. 1. The 15 magnetic vanadium atoms of the 
K6[Vi5As6042(H20)] molecule. They form two hexagonal 
layers and an inner triangular layer of vanadium atoms sand- 
wiched within. The arrows show the lowest energy spin con- 
figuration found from DFT calculations. The six exchange 
parameters used in the Heisenberg Hamiltonian are shown as 
lines. Schematically displayed are energy levels of the Kramer 
doublet (5*= 1/2) ground state and the low lying quadruplet 
(5=3/2) separated by A. 

FM alignment. Thus a FM coupling is obtained if the 
V-0 hopping matrix elements into the higher d-orbital 
are significantly larger than the matrix elements for the 
hopping of O electrons into the lowest energy d-orbital. 
The occurrence of such interactions are possible in a low- 
symmetry system such as the one studied here. 

Even with this FM interaction, our spin Hamiltonian 
yields an 5=1/2 ground state composed largely of Ising 
configurations similar to the one depicted in Fig. 0. This 
Ising configuration was predicted from the J's from ear- 
lier fits to DFT energies and corresponds to the ground 
state DFT configuration (I). 

We have fully diagonalized the Heisenberg Hamilto- 
nian (|^). Using all symmetries, the largest irreducible 
many-body subspace has dimension 2145. We find a spin- 
1/2 Kramer doublet as the ground state with a low-lying 
spin-3/2 quadruplet as shown in Fig. |l|. The rest of the 
spectrum is well separated from these eight states. The 
large value of J binds the spins in the hexagons into sin- 
glets. The low-energy physics arises from the inner tri- 
angle spins interacting with each other both directly and 
with an effective coupling through the hexagons, yield- 
ing the doublet-quadruplet spectrum jl^ . There are two 
important energy scales in the spectrum: A, the gap be- 
tween the doublet and quadruplet, and J, the energy at 
which the singlets in the hexagon break and the molecule 
starts to behave as more than three spins. 

The low-energy effective interaction between inner tri- 
angle spins proceeds via Jl and J2 (which frustrate each 
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other) to a hexagon singlet, via J' and J" to a neigh- 
boring singlet, and finally via Jl and J2. A larger con- 
tribution to A comes from a direct interaction through 
J3 mediated by hopping through O and As levels. Thus 
simple perturbation theory [|| yields 



^_i {J2-Jl)\j"-J') 
4 J2 



(2) 



Comparing our calculated susceptibility with experi- 
ment [ p^ , we find the low-temperature behavior indi- 
cates our doublet-quadruplet gap A « lOK to be signif- 
icantly larger than the experimental value of A w 3.7K, 
while the high-temperature behavior shows our calcu- 
lated value of J is too large. Both of these discrepancies 
can be explained almost entirely by J's that are too large 
within the density- functional-based treatment, as known 
for other vanadium systems 

Agreement with experiment for the low temperature A 
can be achieved by dividing all J's by a constant factor 
of 2.9. Fig. H shows our calculated result with rescaled 
J's compared to experimental data from Chiorescu. A 
uniform scaling of our calculated exchange parameters is 
not able to obtain the right low and high temperature 
behavior at the same time. The high temperature be- 
havior could be improved by futhcr reducing J with cor- 
responding adjustments to the other couplings to keep A 
constant. 

A set of AF interactions |^ also fits the experimental 
results. In fact any set of parameters with the correct 
values of J and A given by Eq. (||) will fit the effective 
moment well. To support future efforts aimed at experi- 
mentally confirming our relative parameters, we have cal- 
culated the spin-spin correlation function Cs = A{Sl-S'j), 
with q representing the arbitrary quantization axis. Since 



TABLE I. DFT energies [E) of calculated Ising configura- 
tions, energies obtained from the fit, and 4(5' SJ) along each 
of the six bonds. Also included is the anisotropy shift 5 for the 
Ms = S state of each Ising configuration. A least square fit of 
this data leads to exchange parameters of J=290.3, J'=-22.7, 
J" = 15.9, Jl=13.8, J2=23.4 and J3=0.55 meV. 
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FIG. 2. The effective moment ^^tt = \/3xT calculated with 
all J's uniformly rescaled by dividing by a constant factor of 
2.9. The diamonds are the actual experimental values from 
Chiorescu. The inset shows the low temperature behavior. 



J is always largest, the spins tend to form a singlet 
along this bond, yielding Cs{J) = —1 (AF correlation). 
The three inner triangle spins form a non-coUinear ar- 
rangement with Cs — —0.33 between all pairs of spins. 
The direct J3 coupling reduces the correlations between 
triangle and hexagon spins. We find C5(J1) = -1-0.02 
and C5(J2) = —0.02 (FM and AF correlations, respec- 
tively), while the coupling constants of Ref. 1^ yield 
Cs(Jl) = +0.12 and Cs(J2) ^ -0.14. This prediction 
should be measurable with neutron scattering. 

Chiorescu et al. observe that rotation of the spin pro- 
jection is achieved without encountering a barrier |l^] 
and Dobrovitski et al. posit that the V15 molecule is 
indeed a low anisotropy system [p^ . As shown below 
the existence of either easy-plane or easy-axis anisotropy 
will shift the Ms=3/2-l/2 Zener-Landau tunneling tran- 
sitions that have been observed by Chiorescu et al. jl^ . 

Recently, Pederson and Khanna have developed a new 
method for accounting for second-order anisotropy ener- 
gies This method relies on a simple albeit exact 
method for spin-orbit coupling and a second-order per- 
turbative treatment of the spin-orbit operator to deter- 
mine the dependence of the total energy on spin pro- 
jection. Initial applications to Mni2 lead to a density- 
functional-based second-order anisotropy energy of 55. 7K 
which is in essential agreement with the experimen- 
tally deduced values of 55. 6K [^. We have generalized 
this methodology to systems with arbitrary symmetry 
and have calculated the anisotropy energy for several dif- 
ferent spin configurations of the V15 molecule. 

We have calculated the anisotropy energy for the low- 
est Ising configurations with one, three, nine and fifteen 
unpaired electrons, as given in Table |. In all cases we 
find that the V15 possesses weak easy-plane anisotropy. 
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This result ensures that anisotropy effects will not change 
the total spin of the V15 ground state. Examination of 
the expression for the second order anisotropy energy 
in Ref. |16|, shows that such energies do not necessar- 
ily scale as the square of the total moment. Indeed, as 
shown in Table |, we find that the anisotropic effects are 
in fact only weakly dependent on the total spin and that 
the energy of the Ms=S state increases by approximately 
^=0.8K to 1.8K. Chiorescu et al. show that the broaden- 
ing of the Zener-Landau tunneling fields decreases with 
temperature for the |l/2, 1/2) to |l/2, —1/2) transitions 
but are independent of temperature for the |3/2, 3/2) to 
the |3/2,l/2) transitions |jl^. This behavior is exactly 
what is expected from a sample containing weak easy- 
plane spin anisotropy. At sufficiently low temperatures, 
only the S = 3/2 and S = 1/2 states are relevant, and 
the field-dependent crossing of these states depends on 
whether the magnetic field is parallel or perpendicular to 
the easy-plane. The broadening is proportional to the dif- 
ference of the magnetic anisotropy energy for the different 
spin configurations involved. Due to the small anisotropy 
in Vi5 the effect will be small. Although it is not pos- 
sible to translate the DFT obtained anisotropy energies 
directly to the quantum mechanical many-spin ground 
state discussed here, we obtain from these energies tun- 
nel field broadenings between 0.1 T to 0.48 T which en- 
velope the experimentally observed field broadening of 
about 0.2 T ||l^. In powdered samples, the small easy- 
plane anisotropy would lead to a broadening in the tun- 
neling field and in single crystals the effect would change 
the tunneling fields as a function of field orientations. 

To summarize, we have performed accurate all-electron 
density-functional calculations on the V15 cluster as a 
function of geometry and spin configuration. By dy- 
namically coupling the mean-field density-functional ap- 
proach to exact-diagonalization of a many-spin Heisen- 
berg representation, we have efficiently determined the 
lowest density-functional configurations and the entire 
Heisenberg spin excitation spectrum. Our calculations 
suggest that the small experimentally observed orienta- 
tional dependence of the tunneling field for the Ms =3/2 
to Ms = 1/2 is a signature of configuration dependent 
magnetic anisotropy in this molecule. The method used 
here is general and allows one to characterize both sys- 
tems which are potentially useful for magnetic storage 
(Mni2, Feg) and systems which show quantum coherence, 
such as the one studied here. 
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